Nanoscale Smoothing and the Analysis of Interfacial Charge and Dipolar Densities 
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£S) , The interface properties of interest in multilayers include interfacial charge densities, dipole den- 

sities, band offsets, and screening-lengths, among others. Most such properties are inaccesible to 
direct measurements, but are key to understanding the physics of the multilayers. They are contained 
within first-principles electronic structure computations but are buried within the vast amount of 
quantitative information those computations generate. Thus far, they have been extracted from the 
numerical data by heuristic nanosmoothing procedures which do not necessarily provide results in- 
dependent of the smoothing process. In the present paper we develop the theory of nanosmoothing, 
establishing procedures for both unpolarized and polarized systems which yield interfacial charge 
^ ■ and dipole densities and band offsets invariant to the details of the smoothing procedures when the 

criteria we have established are met. We show also that dipolar charge densities, i. e. the densities 
of charge transferred across the interface, and screening lengths are not invariant. We illustrate our 
procedure with a toy model in which real, transversely averaged charge densities are replaced by 
sums of Gaussians. 

PACS numbers: 73.30.+y,73.40.Rw,73.40.Qv,73.90.+f 

S: 

^ \ I. INTRODUCTION 

a 



One of the central problems of physics from the mid nineteenth century on has been how to make the transition 
from a microscopic to a macroscopic theory of matter. This problem has two aspects. The first is the task of deriving 
the equations governing the macroscopic behaviour of matter from the underlying microscopic equations, exemplified 
by the derivation of Maxwell's macroscopic equations from the microscopic theory of charges and fields in vacuum. 
This problem is elegantly solved by a coarse-graining procedure which takes an average over "physically infinitesimal" 
regions. Such an elementary region is chosen to be small enough to let the average of a quantity follow all the changes 
| that are observable at the macroscopic level, but large enough compared with characteristic atomic dimensions for 
7-H i it to contain so many particles that the behaviour of an individual particle has a negligible effect on the average 
quantity. This coarse-graining procedure smooths over the atomic-scale fluctuations in physical quantities, leaving 
only the slow spatial variation of their macroscopic components. A beautifully clear derivation of the macroscopic 
Maxwell's equations, first derived by H. A. Lorentz in 1902, 1 can be found in Rosenfeld's Theory of Electrons, a 
now forgotten classic^ In the coarse-graining procedure there are three distance scales: Ai, the scale on which the 
macroscopic quantities vary; A2, the scale on which the smoothing is carried out; and A3, the microscopic scale, that 
is, the atomic scale. For the procedure to work, Ai must be sufficiently larger than A3 that the pair of inequalities 
■ Ai >> A2 >> A3 can both be satisfied. Indeed, this criterion makes clear the distinction between macroscopic and 
' microscopic. 

The second aspect is the task of deriving macroscopic constitutive equations from microscopic properties. An 
O ■ early example is the derivation of the Claussius-Mossotti equation^ which provides the link between the microscopic 
J> , polarizability (response of the atoms or molecules to the local electric field) and the macroscopic dielectric constant. 

The systems of interest in the present paper are those with interfaces between quite different materials. A planar 
r> ] capacitor comprised of an insulating layer sandwiched between metallic electrodes is a good example of such a system. 
Further examples are heterojunctions between different semiconductors and Schottky barriers at semiconductor-metal 
interfaces. In all such systems, there are multiple causes of charge transfer across or to the interfaces. These can 
include charge transfer to establish spatial uniformity of the chemical potential (the Fermi level), charge transfer in 
response to the local change in chemical composition across the interface and to the interface-induced atomic relaxation 
of the structure, charge accumulation to screen the interface charge density associated with the termination of bulk 
polarization at the interface, and charge accumulation attendant to charging or shorting 5,5 of a capacitor. Interface 
dipole densities arise from such transfer of charges and are responsible for offsets of the average electrostatic potentials 
across the interfaces, a dominant factor in determining Schottky barriers and valence and conduction band offsets 
in semiconductor heterojunctionsi 7 ^ In all such cases, the charge density of each material is perturbed, and the 
perturbation is localized near the interfaces, usually within a few interatomic distances. 
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The case of thin ferroelectric films between metallic electrodes provides a good illustrative example of this general 
class of systems. Since the early seventies a phenomenological model has been developed^ ' 10 ' 11 ' 12 to explain the 
modification of the polar phases (substantial reduction of the spontaneous polarization for small thicknesses, or even 
the complete suppression of ferroelectricity below a certain critical thickness) and of their thermodynamic properties 
(depression of the transition temperature with respect to that of the bulk material). The model, mainly due to Batra 
and coworkers, relies on three basic assumptions: (i) the polarization charge lies in a sheet right at the interface, (ii) 
the surface polarization charge density equals the magnitude of the polarization inside the thin film, and (Hi) the 
free compensation charge spreads out at least over a finite distance A within the electrode, decaying exponentially 
towards its interior as in the Thomas-Fermi approximation. In this model, the screening length A is dependent only 
on intrinsic properties of the electrode, such as the density of free carriers or the dielectric constant. All effects are 
neglected which might come from a particular choice of the electrode/ferroelectric interface, such as the different 
chemical bondings formed at the junction or the interpenetration of the electrode and dielectric/ferroelectric wave 
functions that might screen the polarization charge in part within the insulator, reducing therefore the magnitude 
of the interface dipole density. Atomic level charge fluctuations are neglected as well, implying a smoothing on a 
microscopic scale distinct from the smoothing on the coarse-graining scale in the derivation of macroscopic equations 
and properties. 

The need to go beyond such simple models is well illustrated by this interpenetration of wave functions across 
the interface. As pointed out first by Heine^ and later by Tejedor et al^ and Tersoff^ the bulk Bloch-states of 
the metal with energies below the Fermi level of the metal and within the semiconductor band-gap and its valence 
band decay exponentially inside the semiconductor (and, indeed, might have a significant amplitude for a few layers 
from the interface), creating a continuum of gap states [the so-called metal-induced gap states (MIGS)]. Achieving 
deep understanding of interface properties with quantitative predictive power and free of adjustable parameters when 
they are determined by such effects occurring at the atomic scale requires first- principles computations. In recent 
years it has become possible to carry out first-principles calculations for systems of the complexity of those under 
discussion here. These simulations provide a wealth of information at the atomic level about the structural and 
electronic properties of materials and their responses to various external perturbations^ Some quantities, such as 
the microscopic charge density distribution p (f) or the corresponding electrostatic potentials, are routinely available 
from first-principles calculations. 

The question becomes how to extract from the immense detail provided by the first-principles computations reliable 
values of the physical quantities of interest — interface charge and dipole densities, screening lengths, etc. — which 
enter the pseudo-macroscopic models currently used. Two major difficulties arise. First, coarse graining is inapplicable 
because the relevant distance scale for interface properties is the atomic scale, i. e. A2 ~ A3. Second, the charge- 
density changes associated with maintaining the constancy of the Fermi levels can be orders of magnitude smaller than 
the unperturbed bulk charge densities, themselves very rapidly varying functions of position, reflecting the underlying 
atomic structure. The relative magnitudes of these changes in the microscopic charge densities are illustrated in Fig. [T] 
Moreover, the polarization-induced charge densities and their screening charge densities can be smaller by additional 
orders of magnitude than those arising from imposing Fermi-level constancy. Therefore, all the interface-related dipole 
densities are overwhelmed by much larger variations of the total microscopic charge density. 

It is still possible to carry out smoothing of the computed charge density at the nanoscale where the condi- 
tions for coarse-graining are not met, providing sufficient care is taken. A heuristic smoothing procedure has been 
introduced ^ 17 ' 18 to extract the quantities of interest from the results of first-principles charge-density computations 
as described in more detail below. However, that procedure has not yet been systematically analyzed to establish 
the conditions under which it accurately extracts the quantities of interest: surface charge densities, surface dipole 
densities, etc. In the present paper, we introduce such an analysis. In addition, in order to focus only on the per- 
turbations introduced by the nanosmoothing procedure, avoiding other sources of numerical noise coming from the 
first-principles simulations, we illustrate the analysis with a toy model whose accuracy can be arbitrarily improved. 
The application of this theory to first-principles computations on realistic ferroelectric capacitors is the subject of a 
forthcoming paper. Despite this focus on the ferroelectric capacitor, our analysis is of general utility for the extraction 
of interface properties for all multilayer systems. 

The rest of the paper is organized as follows. We set the grounds of our discussion of simulations of interfaces 
from first-principles and define the microscopic behaviour of the different quantities that are the targets of our study 
in Sec. HU In Sec. IIIII we define the interface quantities of interest for an unpolarized interface. In Sec. IIVI 
we describe the difficulties encountered in defining precisely the location of a reference interface. We develope the 
theory of nanosmoothing of unpolarized systems in Sec. |V]with particular attention to questions of the sensitivity of 
quantities of interest to the smoothing procedure. We describe in Sec. I VII a toy model used to illustrate all of the 
essential elements of the smoothing theory. In Sec. I VIII we present the results of the toy model computations for the 
unpolarized case. We generalize the nanosmoothing theory of Sec. |V]for polarized systems in Sec. I Villi and present 
the corresponding toy model results in Sec. [TX] Finally in Sec. [}Q we summarize our results, emphasizing the specific 
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FIG. 1: Charge density laterally averaged parallel to (001) planes [see Eq. ([1} below] of bulk BaTiOa (a), SrRuOs (b), and 
(c) an unpolarized planar capacitor comprised of 4 layers of BaTiOs between metallic electrodes made of 5 layers of SrRuOs, 
all from first-principles. The charge-density profile of the capacitor looks like a juxtaposition of the bulk charge densities of 
the two materials, highlighting the fact that the charge transferred from one material to the other to establish the constancy 
of the chemical potential is overwhelmed by much larger variations of the total microscopic charge density. Details of the 
first-principles calculations can be found in Refs. andH|. The unit of the charge density is electrons/bohr 3 . 



criteria smoothing functions must satisfy to yield interfacial properties insensitive to their parameters and specifying 
which of our results are new. 



II. SIMULATION OF INTERFACES FROM FIRST PRINCIPLES. 



First-principles calculations of interfaces between two materials, where there is no periodicity in at least one direc- 
tion, are almost universally done by means of the supercell approximation. 19 Within this approach a basic unit cell 
that contains a suitable number of multiatom layers of the two materials is periodically repeated over all space (Fig. 
[2]). For the interfaces within a nanostructured multilayer to be well defined, with properties distinct from those of 
the bulk-like regions between them, the widths of the layers of each material introduced in the construction of the 
basic unit cell must be large enough to avoid the interaction between adjacent interfaces through the bulk materials, 
so that the calculation accurately represents an isolated interface. 

Throughout this work, we shall assume that the interface is oriented along the z axis, and each material is periodic 
in the plane parallel to the interface, referred to as the (x, y) plane. 

Precisely the same methodology is used to treat nanostructured multilayer materials, and in this paper we do not 
distinguish between the cases, specializing to multilayers in which the individual material thicknesses are large enough 
for the interfaces to be noninteracting. 

The microscopic charge densities p (r) [see Fig. HJc)] and electrostatic potentials V (r) provided by the first-principles 
computations for the previously described supercells are continuous functions periodically repeated in space with the 
periodicity of the supercell (a sc in Fig. [5]). A few interatomic distances away from the interfaces, the microscopic 
quantities recover their bulk features. In other words, if the layer widths are large enough, it is possible to identify 
"bulk-like" regions in the middle of each of the layers that constitute the superlattice, with large variations of the 
microscopic quantities with the same periodicity as in the bulk, reflecting the underlying atomic structure. 

Many of the interface-related quantities that we shall define below depend only on what happens in the direction 
perpendicular to the interface, where the actual discontinuities of the physical structure occur. Since we assume in- 
plane periodicity, we can trivially eliminate the in-plane dependence by taking a planar average of the corresponding 
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FIG. 2: Schematic view of a supercell used in first-principles simulations of interfaces. The z axis is taken as normal to the 
interface, whereas the plane parallel to the discontinuity is taken as the (x,y) plane, m and ri2 multiatomic layers of two 
materials, whose bulk lattice constants in the out-of-plane direction are respectively a\ and 02, are stacked to build a basic unit 
cell (hatched rectangle) that is periodically repeated in space. a sc is the length of the supercell in the out-of-plane direction, 
with a sc — max + 11202. In the figure, n\ — 2, and 112 = 3. 



microscopic quantity, e. g., for the electron density, 



P ^ = s J J s p ^ dxd y-> 



(1) 



where S is the area of the interface unit cell. 

As is proven in Appendix ^ this transverse averaging has no effect on the Poisson equation, which reads after 
averaging 



V 2 V (r) = V 2 V (z) 



d 2 V(z) 
dz 2 



-Airp (z) . 



(2) 



III. INTERFACE QUANTITIES OF INTEREST; THE UNPOLARIZED CASE. 



One of the most important physical properties of heterojunction devices is the band offset or Schottky barrier at an 
interface, that is, the relative positions of the energy levels on both sides of the interface. In the case of semiconductor 
heterojunctions, the valence-band offset (VBO) [conduction-band offset (CBO)] is defined as the difference between 
the positions of the tops of the valence bands [the bottoms of the conduction bands] of the two materials. In the case 
of metal-semiconductor contacts, we can define the p-type (n-type) Schottky barrier as the difference between the 
Fermi level of the metal and the top of the valence band (bottom of the conduction band) of the semiconductor. These 
differences in the band positions determine the effective barrier for electron or hole transport across the junction. 

The computation of such effects from first-principles cannot be achieved by a direct comparison of the corresponding 
single-particle energies (tops of the valence bands, bottoms of the conduction bands and/or Fermi level of the metal) 
in the two compounds as obtained from two independent bulk band-structure calculations. The reason is the lack of an 
intrinsic energy origin to which to refer all the energies: in a first-principles simulation, the Hamiltonian eigenvalues 
are referred to an average of the electrostatic potential that is ill defined for an infinite system- - where, due to 
the long-range nature of the Coulomb interaction, it is defined only to within an arbitrary constant. Consequently, 
together with the eigenvalue difference, we must consider both the shift of this average between the two materials and 
the redefinition of the averaging process so that it is appropiate to the multilayer system under consideration. As was 
mentioned in Sec. [H the coarse-graining procedure used historically fails when applied to the results of first-principles 
computations. We designate the averaged potential as (V^) and its shift as A(V U ). The brackets () indicate that the 
averaging process is not yet defined. We shall label all quantities of physical interest of the unpolarized system by a 
subscript u. This potential shift depends on the dipole induced by the electronic charge transferred from one side of 
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the interface to the other after interfacial hybridization. As the charge transfer depends not only on the materials that 
constitute the interface, but also on intrinsic interface effects such as the chemical composition (termination of each 
material at the interface), on particular orientation and on other structural details, the shift can only be obtained 
from a self-consistent calculation on a supercell including both materials. This ensures that the averaged electrostatic 
potentials of both materials in the "bulk-like" regions defined in the previous section, where all the physical quantities 
recover the bulk features, are expressed with respect the same reference and allows a direct extraction of the shift. 

This shift of the averaged electrostatic potential should be directly related to an averaged interface dipole-moment 
density (p„), 



A(K) = 4tt(p u ). (3) 

We prove in Appendix [Cl that Eq. ([3|) holds for the specific definition of the average procedure () introduced in Sec. 
IV Al below and discussed in this context in Sec. IV Bl 

From a fundamental point of view^i the charge Q u , and electric dipole moment at the interface p u are defined 
respectively as the zero and the first moment of the total microscopic charge density, 



Q u = I dz p u (z) , (4a) 

J zi 

dz zp u (z) . (4b) 

-i 

For the previous equations to be meaningful and truly represent an interface quantity, the thicknesses of the adjacent 
layers must be wide enough so that they contain regions within which the charge density is essentially unaffected by 
the presence of interfaces within which Z\ and zi could be located. 

However, from a practical point of view, such a definition poses serious questions. Indeed, both Q u and p u are 
ill-defined due to the large and rapid oscillations of the microscopic charge density [see Fig. (He)]. In fact, different 
choices of the integration limits z\ and zi yield widely different values of the charge and interface dipole moment. 
Only in the extreme case of a Clausius-Mossotti model, in which the total charge is unambiguously decomposed 
into an assembly of localized and neutral charge distributions, so that a unit cell can be chosen with no charge at 
the surface, the dipole moment of a periodic charge distribution would be well-defined as the integral of the first 
moment of the charge density. However, any Claussius-Mossotti approach does not correspond to reality, particularly 
in materials where delocalized covalent charge is present^ This ambiguity in the definition of p u with respect to the 
boundaries of the region within which the dipole moment is computed is closely connected to the problem of defining 
the polarization of a periodic system from the charge density^ 



IV. THE DIFFICULTY OF DEFINING A REFERENCE INTERFACE. 

In order to get rid of bulk effects and extract interface-related features, some author a 24 > 25 ' 26 i 27 have defined an ideal 
interface by stacking alternate slabs, each of them made from slicing the planar average of the bulk charge density of 
the corresponding material perpendicular to a particular direction. Let us define pp 1 ' (z) and (z) as the planar- 
averaged [Eq. ([I])] bulk charge densities of the left and right material respectively, unaffected by the presence of the 
interface. They are locally periodic in z, and, assuming that both bulk materials are unpolarized, the net charge and 
the dipole moment vanish within each bulk unit cell, 



M 

bu.lk 



dz -p[ s) (z) = 0, (5a) 



w dz zp^ (z) = 0, (5b) 

bulk 

(s) 

where a bulk is the lattice constant of the bulk unit cell of each material in the z direction, and s refers to the side of 
the interface considered, 1 or 2. The bulk unit cell boundaries are chosen so as to preserve inversion symmetry. 

Replicating the bulk charge densities up to the as yet unspecified interface from each side, we could define the 
reference charge density p for all z as 
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Po ( z ) = { 2 



Po ^ ( z int) + Po ^ (Zint) 



(6) 



where z,- n < is the coordinate assigned to the position of the interface. p (z) is discontinous at the interfacial plane 
Zint by its very definition. 

We now define the interface-induced deformation of the charge density Ap u (z) as 



A Pu 0) = Pu ( z ) - Pa ( z ) ■ ( 7 ) 

If the thicknesses of the two layers are wide enough, A~p u (z) becomes negligibly small over the ranges R\ in the left 
material and R 2 in the right material [see Fig. UHb)] . The interface region can thus be identified as comprised of those 
ranges where A~p u (z) differs significantly from zero, in other words where the microscopic charge density differs from 
the relevant bulk values. The interface charge and dipole density associated with A~p u (z) are defined as, 

AQ U = dz Ap u (z) , (8a) 

Jzi 

A Pu = dz zAp u (z) . (8b) 



The advantage of this approach is that both AQ U and Ap u are well defined quantities with respect the location of 
the integration limits z\ and z 2 in Eqs. ((8a|) - (f8b|) . provided that Z\ lies in R\ and z 2 lies in R 2 . 

However, this approach has pitfalls. In particular (i) the position of the interface, Zi n t m Eq. ©, is not yet 
specified. Some recipes have been given for how to cut the bulk slabs, but they have limited applicability. One case 
is for common anion heterostructures with non-relaxed interfaces along high-symmetry planes, such as the (001)24, 
(110) 25 , or (111)2£ interfaces of GaAs/AlAs superlattices. As soon as an interface-induced rippling of the atomic 
layers is introduced, for instance after an atomic relaxation of the interface geometry, the problem of defining the 
position of the interface worsens, (ii) Therefore, the interface charge and dipole densities are not unique, since they 
depend critically on where the mathematical surface representing the interface is chosen, and no objective criterion 
for locating it has been established. In particular AQ U and Ap u [Eqs. (|5a)) - (|8b[) ] equal Q u and p u [Eqs. (T4"a)) - (l4b[) ] 
if and only if Zi nt — z\ contains an integer number of unit cells of the left material and z 2 — Zi nt contains an integer 
number of unit cells of the right material. The location of the interface determines where on one side, the bulk charge 
density of the left material is subtracted, and on the other side that of the right material. A different choice of 
the mathematical interface can produce very different charge and dipole densities. In addition, comparison between 
different interface orientations makes little sense with this definition. 7 (Hi) Since the interface dipole moment is 
dependent on the reference charge density, the corresponding potential drop at the interface (Arfj po ; e in the notation 
of Refs. [24.25.26.27) must be too. However, it is important to note that the potential drop generated by Ap u (z) is 
only part of the total potential shift. The total charge density of the interface is given by 



Pu 0) = Po ( z ) + A Pu ( z ) » (9) 

thus the potential shift associated with 7j (z) must be included as well. Although the existence of a potential shift 
generated by Po ( z ) is general, we shall explain its origin only in the particular case where Zi n t — Z\ contains an integer 
number of unit cells of the left material, and z 2 — Zi nt contains an integer number of unit cells of the right material. In 
this particular situation both slabs used to construct the reference charge density in Eq. ^ have neither a charge nor 
a dipole moment. Under these circumstances, the potential shift is the difference in the locally averaged potentials 
produced by the p of each material in the regions R\ and R 2 respectively. For this shift to be an interface property, 

the layer width of each material must be large enough that the local average (Vj^ b ) has approached the unperturbed 
value within the center of an slab. The planar averaged potential within any point of the slab will be given by 2 ^ 
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V^{z) = 2ir I dz'\z~z'\pi s \z'), (10) 

so the local average can be computed as 



(O = -S- / dz [dz'\z- z'\fc\z). (11) 

®bulk central cell J 

Since (V^ b ) depends on the charge density distribution, it differs for the left and right slabs used in the construction 
of the reference charge density and produces an additional shift. Therefore, the total potential drop at the interface 
is the sum of the potential drop generated by Ap u (z), Adi po ie, and the difference of the average potential of the two 
reference slabs A re /. Only the sum A re f + Adi po ie is independent of the reference charge density chosen and is a 
physically measurable property of the interface. Since each term in the sum is sensitive to the arbitrary location of 
the interface, each must be computed accurately enough for the sensitivity to disappear from the sum. As the charge 
density shifts of interest are so small, this is an unnecessary burden, removed by the use of a proper nanosmoothing 17 
procedure as shown in Sec. [V] 



V. NANOSMOOTHING. 



A. The procedure. 



A procedure to eliminate charge fluctuations in the regions of the material which do not contribute to the interfacial 
hybridization, thereby localizing the physically relevant charge densities to the interface, consists of filtering out the 
periodic oscillations of microscopic quantities, which typically follow the underlying atomic structure, preserving only 
those features that emerge in the vicinity of the interface. 

To obtain this smoothed charge density, we have followed the recipe given by Baldereschi et al. in Ref. [I?] and 
generalized by Colombo and coworkers for lattice mismatched heterostructures in Ref. Il8l . Starting from the planar- 
averaged charge density ~p u (z), we construct the smoothed density ~p u (z) by convoluting it with a smoothing function 

/(*), 



Pu i z ) = I dz f(z- z )p u {z ) (12) 

which has the following properties: 28 



/(*)>0, 
f(z) = 0, 
f(-z) = .f(z), 

dzf(z) = 1, 



\z\<L, 
\z\>L, 



(13a) 
(13b) 
(13c) 

(13d) 



so that 



Pu 0) 



r z+L 



dz f(z - z )p u {z ). 



t-L 



(14) 



In addition, / (z) should be monotonic in \z\ and sufficiently smooth itself. L should be chosen on the scale of the unit 
cell length or larger, but smaller than the widths of the left and right layers. A sharper criterion for L is introduced 
below. We define this as the averaging procedure () left unspecified above in Sec. Mil 

The particular smoothing function we have used, following Refs. |$ll2landll8|, is the convolution of two square- wave 
filter functions: 
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f(z- z) = J dz u h (z- z )uji 3 (z -z), (15) 

where 

cJt(z) = le(L-\ z \\, (16a) 
f 1, z > 

0(z) = < (16b) 
[ 0, z<0. 

Giustino and coworker s 29 ! 30 propose convolution with a Gaussian kernel that can be an approximation to the asymp- 
totic limit of a convolution of a large number of square- wave filter functions. This method is best suited for superlattices 
where crystal deviates from perfect periodicity far away from the interface so that it is not possible to define regions 
where the interface-induced charge density vanishes, or in disordered three-dimensional systems with short-range 
order. Even more general functions can be used, providing the criteria established above are met. 
The explicit dependence of / (z), defined in Eq. (fTS")) , on z is: 



|*| > 



h+l-i 
2 ' 



/(*) 



hh [ 



2 

J_ 



(17) 



< \z\ < 



\h-h 




FIG. 3: Filter functions used for the smoothing of the charge density and potential in the present work. Two square wave 
functions of lengths h and h, as denned in Eqs. (|16ap and (|16b|l . are represented in panel (a). The corresponding convolution, 
given by Eq. (|15[) . is shown in panel (b). For numerical calculations, the square- wave functions are expanded in Fourier series. 
The small oscillations in the vicinity of each discontinuity are due to the Gibbs phenomenon^ 



as shown in Fig. [3] Z> is the greater of l\ and fa. 

After nanosmoothing the planar averaged charge density 7j„ (z), the resulting charge density ~p u (z) is a continuous 
function that joins smoothly at the interface. Even though we have not subtracted anv reference charge, the smoothed 
charge density becomes negligibly small over ranges 7?j in the left material and R 2 in the right material [see Fig. 
E{c)], and the interface region can be unambiguosuly defined as the region where the smoothed charge significantly 
differs from zero. Note that R s and R s will differ in general, as shown in Fig. [5] and discussed further in Sees. IV CI 
andEIIl 



9 



B. The Poisson equation and potential shifts. 

Providing the filter function f(z) satisfies the following additional conditions, 

d 2 f(z) 



dz 2 
df(z) 



dz 



exists, 



z=-L 



df(z) 



dz 



= 0, 



z=+L 



/WU L = /WU +L = o. 

The Poisson equation remains invariant after the nanosmoothing (see Appendix [B| and transforms into 



(18a) 
(18b) 
(18c) 



W(z) 



d 2 V(z) 
dz 2 



-4arp (z) . 



(19) 



The nanosmoothing function f(z) defined in Eq. (|16ap and Eq. (|16b[) violates condition Eq. (|18a[) at its end points 
2 = ±|, where f(z) is discontinuous so that condition (|18c[) cannot be unambiguously applied. Similarly, the f(z) of 
Eq. (fl7| violates condition (jTHaJ) at the four points z — ± l1 ^ 1 ' 2 , and df /dz is discontinuous at the points z = ± ll ^ 12 
so that condition (|18b|) cannot be unambiguously applied. Nevertheless, one can think of the f(z) as a distribution, 
a family of smooth functions all of which meet conditions (|18a[) - (|18cj) and which approach the f(z) of Eq. (JT7J) as 
their limit. In practice, because the smoothing operation is a convolution, Eq. (|12p . one carries out smoothing via 
fast Fourier transformations. The family of the finite Fourier series involved is thus a distribution which converges to 
f(z) in the limit, meeting conditions (|18a|) - (|18c|) along the way. 

Eq. (|19p holds in general for both the unpolarized and the polarized cases, so subscripts have been omitted. 
Otherwise the subscript u is used because we are presently treating the unpolarized case. 

We prove in Appendix [Cl that the full electrostatic potential shift AV U is given by the nanosmoothed dipole density 



where 



AV U = 4irp u , 



(20) 



Pu = I dz z Pu ( z ) 
lz 1 



(21) 



C. Insensitivity of the dipole- moment density and potential shift to the smoothing function. 

The microscopic charge density ~p u must return to the bulk microscopic charge densities ~p^ and J>^ in the regions 
Ri and R 2 for it to be possible to ascribe physical properties specifically to individual interfaces. However rapidly it 
approaches those values, the approach must be complete within R\ and i?2, as summarized by Eq. 



fp 1} (z), zeR u 



Pu ( Z ) = < 



g(z)^p^(z), ztRi, 
g(z)^pW(z), z^R 2 , 
pf\z), z&R 2 , 



(22) 



with g (z) such that p u (z) and all its derivatives are continuous (the cusps at the nuclei are washed out by lateral 
averaging). Consequently, locally within R\ and R 2 , the microscopic charge density ~p u can be represented by the 

Fourier transforms of ~p^ and ~p^ , 
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Pa 



K)f> = -pr-; n S %; z£R s ; s = 1, 2. (23) 

"huZfc 

In Eq. ([23)1 Aq vanishes because of charge neutrality. 

The smoothing function f(z — z ), Eq. (|15[) . is a convolution of two square wave functions loi s {z — z ), Eqs. 
(116a|) - (|16b|) . s = 1, 2. The order of the within the convolution is immaterial, so the oji s can be applied to the 
nanosmoothing of 7j first, in the two-step nanosmoothing process implied by use of f(z — z ) in Eq. (|12p . As long as 
l s is an integer multiple of the lattice constant of material s, all contributions to p Q from A n 0, are smoothed to 

zero, leaving only which itself vanishes. However the regions R s within which p u vanishes lie within R s because 
smoothing ~p u within R s brings into p u values of ~p u for z outside R s . The limits z s for the determination of p in Eq. 
(f2"Tj) must lie within R s , and L must be significantly smaller than the width of R s . 

Similarly the electrostatic potential V s (z) can be expressed within R s as a comparable Fourier series, 

V s (z) = Y / B { n s) e^ >z , (24) 

n 

with the Bn^ fixed by the Poisson equation, Eq. J5J), 

B^ = ^,n^0, (25) 

( s ) 

(s) . . ' 

except for Bq which is influenced by the charge density outside of R s . Upon nanosmoothing, all contributions to 
V u (z) for z within R s vanish except that for n = 0, 

V.(z)=BP t (26) 

which is invariant to the smoothing process. 

The potential shift AV U , defined in Eq. (|C2[) as 

AF„ = V u (z 2 ) - V u (zi) , (27) 

is thus invariant to the smoothing procedure, 

AF U =4 2) -B 1) . (28) 
Moreover, according to Eq. (|20[) . the nanosmoothed dipole-moment density p u is invariant as well, 

(29) 



1 „= 1 

47T 4-7T 



? (2) R (l) 



-^0 - B 



D. The transferred charge density and the dipolar density. 

It is of considerable physical interest to establish the value of the charge transferred across the interface, a difficult 
task. A criterion for establishing the position of the interface is needed, the difficulty of which is discussed in Sec. 
IIVI The rapid, large oscillations of ~p u (z) and the relative smallness of the pertinent component of ~p u (z) make using 
it impractical. On the other hand, if one uses a criterion based on ~p u (z), Zi„t can be sensitive to the smoothing 
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function. Nevertheless, we shall attack the problem using ~p u (z) and attempt to overcome the resulting sensitivity to 
the smoothing function of the position of Zi n t and the amount of charge transferred. 
We start by defining two cumulative charge densities 



Q- 0) = dz p u (z ), 
Q+ 0) = / dz J5 u (z ). 

J z 

For the unpolarized case now under consideration, Q u = 0. Thus, as 



(30a) 
(30b) 



l u = Q-{z) + Q + {z), ze(z 1 ,z 2 ), 

L (z) = -% (z) 



(31) 
(32) 



holds for Vz E (zi, Z2), and it is sufficient to consider either one or the other. Define q as the magnitude of the charge 
transferred per unit area of the interface, the transferred charge density. We estimate q as 



q = sup Q ± (z) , 



(33) 



and estimate Zi nt as 



Zint = arg sup Q±(z) 



(34) 



Now both q and Zi n t are sensitive to the choices of l\ and I2 in the smoothing function /. As l\ and I2 increase, 
Lai 1 can reach across the interface from material 1 to material 2 bringing contributions from ~p{z ), z within material 
2, to ~p(z), z within material 1, and vice versa, thus returning part of the transferred charge back to its origin and 



'(z) 



could 



reducing the value of q. Moreover, since p(z) contains components which oscillate strongly with 2, 

develope multiple suprema or maxima as h and I2 increase in multiples of the lattice constants of materials 1 and 2, 
respectively. This would vitiate the utility of the definitions (|3"3"|) and (|3"4")l of the transferred charge q and the interface 
location Zi n t, respectively, should it happen. We have found that it does happen in the toy model described in Sec. 
|VJ and studied in Sec. IVII1 cf. Fig. [8] below, in the case where the interatomic distance remains unchanged across 
the entire superlattice. Accordingly, as a precaution, the smallest acceptable values of l\ and l 2 should be used for 
/ (z), a single lattice constant of each material, an important additional condition on l\ and l 2 - 

If the transferred charge density were concentrated equally on two surfaces at either side of the interface, separated 
by a distance A, a dipole moment density of magnitude qX would be created. Setting qX equal to the actual dipolar 
density p u allows us to define X u as the dipolar length 



Ai< — 



Pu 



(35) 



where we have restored the subscript u to q as we are dealing with the unpolarized case. 



E. Loss of invariance of physical magnitudes of interest with nanosmoothing. 

Now we can ask whether the physical magnitudes of interest, such as the interfacial charge [Eq. (|4a[) ] or dipole 
moment densities [Eq. (|4bj) ] remain unchanged if the microscopic charge density is replaced with the nanosmoothed 
charge density. In other words, if we define the interfacial charge Q u and dipolar densities p u computed from ~p u (z) 

as 

dz "p u (z) , (36a) 

zi 

22 _ 

dz zp u (z) , (36b) 



12 



then the question is whether Q u — Q u and p u — p u for given integration limits Z\ and z-i . 
Replacing Eq. HU) into Eq. (|55a|) . 



rZ2 rZ2 pz-\-L 

Qu= dz K( z ) = dz / dz'f(z- z')p u (z). 

J z\ J Zi J z — L 



(37) 




FIG. 4: (a) The range of integration in the z, z plane according to Eq. (|37[1 . (b) The range of integration according to Eq. 
(|39l) . The direction of the shading lines, horizontal in (a) and vertical in (b), indicates the variable of which the first integral 
occurs, over z in Eq. (|37|l . and z in Eq. (|39|) respectively. 



The region of integration in Eq. ([37]) within the z, z plane is shaded in Fig. [4|a). Then, the integral of Eq. (|37|) can 
be decomposed into the integral on the central square plus the integrals on the cross hatched triangles 



On the other hand, Q is defined as 



22 rZ2 f 

dz I dz f(z - z )p u {z )+ 

Zi Jzi 
z 2 rz+L i i 

dz I dz f(z — z )~p u {z )- 

Z 2 —L J Z 2 

zi+L i-z\ i i 

dz \ dz f(z-z)p u (z). 

zi J z — L 



(38) 



Q 



dz Pu( z ')= / dz Pu( z ) I dzf(z-z), 

z\ J z\ J z —L 



(39) 



where the region of integration in the z, z plane is now shaded in Fig. Hfb). Decomposing the domain of integration 
as for Q u above gives 
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rz 2 


rz+L 


/ dz 




' Z2 — L 




rzi + L 




/ dz 




'Zl 


Jz-L 



pZ2 rZ2 t 

Q u = dz dz f(z- z )p u (z ) + 

pZ2 + L pZ2 

dz I dz f(z — z )~p u (z )+ 

/ Z2 J z — L 

dz [ + dz' f(z-z)p u (z). (40) 

' zi — L J zi 

Therefore, the equality Q u = Q u is verified if and only if the integrals over the shaded regions on either side of the 
line z = z are the same, that is if the following pair of equations hold 

/ / / fZ2 + L r-Z2 t i / 

dz f(z-z)p u (z') = dz / dz' f(z- z')p u (z), (41a) 

J Z2 J Z — L 

dz f(z — z )~p u {z ) = I dz j dz f(z — z )~p u (z ). (41b) 

J zi — L J zi 

If in Eq. (|41a|) we apply the following change of variables z — Z2 — u and z = Z2 + v in the left hand side and 
z = Z2 + u and z = Z2 — v in the right hand side, then Eq. (|41all transforms into 

L rL — u t>L rL — u 

du I dv f(—u — v) ~p u (z 2 +v)= du dv f(u + v) p u (z 2 — v). (42) 

o Jo Jo Jo 

Due to the parity conditions of the filter function, we know that /(— u — v) = f(u + v). But, ~p u (z2 + v) = ~p u (z2 — v) 
if and only if ~p is even about Z2- Completely analogous reasoning can be applied to Eq. (|41bl) and is omitted here. 

Thus, for Q u to equal Q u , it must be possible to find a z\ and a Z2 about which p(z) is symmetric for \z — z\^\ < L. 

We now show for the unpolarized case that this symmetry condition can be satisfied. For a multilattice consisting 
on alternating layers of two different materials s, s = 1 or 2, there are two distinct interfaces within the supercell 

bounded by z\ and z\ + a sc , interface 1,2 and interface 2,1. We define total charge Q sc and Q sc which are the sums 
of the charges associated with each individual interface, 

Qsc = Q12 + Q21 = / dz' p u (z ), (43a) 

J zi 

pzi+a 3c 

Q sc = Qi2 + Q2i = dz %{z). (43b) 

J zi 

Proceeding in analogy with Eq. (|3"5|) and Eq. (|4"0")) . we obtain 

Q ac = dz / dz f(z-z')p u (z')+ 

J Zi J Zi 

rZi-\-a sc pz-\-L 

dz / dz f(z — z )p u (z )+ 

f zi-\-a sc — L J z\-\-a sc 

dz / dz' f(z-z)p u {z), (44) 

'21 •/ z — L 

Qsc = I dz dz f(z- z )p u (z )+ 

J Zi J Zi 

rZi-\-a sc -\-L pZi-\-a sc 

dz / dz f(z - z )p u (z )+ 

' zi-\-a sc J z — L 

rzi rz-\-L f f 

dz / dz f(z - z )p u (z ). (45) 

' zi —L J zi 
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The fact that p u (z) is periodic in z with period a sc allows as to rewrite Eqs. (|44| and (I45|) so as to establish the 
equality of Q sc and Q sc , 

Qsc= dz / dz f( z - z )Pu( z ) + 

J Z\ J Z\ 

pzi rz+L i i 

/ dz / dz f(z - z)p u (z) + 

Jzi—L J Z\ 

[ 1+ dz [ 1 dz' f(z-z')p u (z') = Q sc (46) 

J Z\ J z — L 

Since the supercell is electrically neutral, so must the nanosmoothed supercell be, 

Q sc = = ^ sc . (47) 

Consequently, from Eq. (|43b[) it follows that 

3l2 = -^21, (48) 

implying that there would be a smooth electrostatic field within each layer if Q 12 is nonzero. The existence of such a 
field would polarize the system in contradiction to the initial condition that the system is unpolarized. We conclude 
that 

^12 = ^21 = 0. (49) 

Thus, for the interface charge to be invariant to nanosmoothing, that is, for Q 12 = Qvi to hold, z\ and z 2 must be 
positioned in R\ and R 2 so that Qi 2 vanishes in the unpolarized case. To do this, one could make an arbitrary choice 
of z\ in R±, say, and then integrate ~p u (z) from z\ up to some z 2 in R 2 at which the integral vanishes. There is no 
need to do this, as it is the nanosmoothing quantities themselves which are of interest. 

Repeating the reasoning for the dipole moment density, we arrive at the conclusion that for p u = p u , the following 
condition must be satisfied 



C Z2 f-z+L f f i pz 2 + L rz 2 r r 

I dz z dz f{z — z )~p u {z ) = I dz j dz z f(z — z )~p u (z ), (50) 

J Z2 — L J Z2 J Z2 J Z — L 

Applying the same change of variables as before, that is z = z 2 — u and z = z 2 + v in the left hand side and z = z 2 + u 
and z = Z2 — v in the right hand side, then Eq. (|50p transforms into 

/>L />L — u />L />L—u 

/ du (zz — u) / dv f(-u-v) p u (z2+v) = du dv f (u + v) p u (z 2 — v) (z 2 -v). (51) 

Jo Jo Jo Jo 

Even in the case of a function ~p u symmetric around z\ and z-2 , the previous condition does not hold in general, and 
the difference between p u and p u amounts to 

/>L r-L — U 

du dv (u — v) f(u + v) p u (z 2 +v), (52) 

Jo Jo 

plus a similar term that comes from the difference in the lower triangles in Fig. BJa) and Fig. [4jb). This has to be 
evaluated for each particular case. 

In conclusion, the interface dipole-moment density is not invariant to nanosmoothing and the interface charge density 
can be made so only by exquisite case in the choice of z\ and z 2 . This is of no concern, as it is the nanosmoothed 
quantites which are of physical interest. 



15 



VI. DESCRIPTION OF THE TOY MODEL. 



As highlighted in the introduction, the components of the density which give rise to the interface-related dipolc 
densities are nearly obscured by the much larger variations of the total microscopic charge density. This atomic- 
scale charge density, routinely provided by any density-functional-based first-principles code, is affected by numerical 
noise and convergence problems inherent in some of the standard approximations in the practical implementations 
of density functional theory (DFT). Therefore, the accuracy of the computations required for extracting the actual 
charge transferred from one side of the interface to the other must be high enough so that the numerical noise of the 
calculations is orders of magnitude smaller than the relevant interface-related charge densities. In this paper, in order 
to illustrate all of the essential elements of the theory of smoothing while avoiding these practical problems, we shall 
define a toy interface model that resembles closely a realistic multilayer material but whose computational accuracy 
can be systematically improved. 

The requirements that such a toy model should meet are: (i) its electron density must be a continuous function; 
(ii) far away from the interfaces, where the interface-induced perturbation of the charge density becomes negligible, 
the toy electron density must tend to two distinct periodic functions on the left and on the right of each interface, 
mimicking the differing behaviour at the bulk level of the materials that constitute the multilayer system containing 
the interfaces; and (iii) the interlayer spacing at the interfaces should be distorted with respect those at bulk so as to 
simulate better the interface induced relaxations that happen in real interfaces. 

In the toy model we propose here, we represent only the laterally averaged density, a one-dimensional function. As 
before, the direction perpendicular to the interface is referred to as the z axis. 

We shall consider atomic-like charge densities gf (again as before, the superindex s = {1,2} refers to the side of 
the interface, left or right, where a given "atom" i is located) of the sum of two gaussians, centered at each atomic 

site zf , 




The positive (negative) gaussian, whose standard deviation is denoted by a n (cr e ), mimics the nuclear charge density 
(electronic charge density) of atom i centered on position Zi. By imposing a n < a e we make the "nuclear" charge 
more confined than the "electronic" charge. Since the gaussians are normalized, a net charge per atomic site can be 
simulated by making A™' 6 ^ A^ s . The parameter 8^ allows us to displace the electronic clouds with respect to the 
nuclei and thereby produce a net dipole moment on a particular atom. 

The atomic-like charge densities are arranged in bulk unit cells that might contain one single atom or a more 
complicated polyatomic basis. In this Section we shall assume that the bulk unit cells of each of the materials that 
form the superlattice contains a single atom that does not carry any charge (A™' s = A^' s V i) or dipole moment 

(<jj — V i). For the polarized case, we refer the reader to Sec. HXl Thus, for the unpolarized interface each material 
can be considered as a one-dimensional monoatomic chain, where consecutive "atoms" are separated by a distance 

cS s \ The interatomic distance at the interface, o m i, is taken as a,i n t = - — — • Then, consecutive interatomic 
distances from the interface evolve smoothly towards the bulk value as a function of the distance to the interface. In 
our simulations, when we move from the interface towads material s,s — {1,2}, the second interatomic distance is 
set up to + ja,( s \ where is the lattice constant of the other material. The bulk value is recovered only at 

the third interatomic distance from the surface. 

A supercell is then built, as described in Sec. |TTJ The basic unit cell, periodically repeated in space, contains a 
suitable number N± and N2 of bulk unit cells of the two materials. The microscopic charge density, ~p(z), is defined 
as the superposition of all the atomic-like charge densities 

Ni N 2 

pW = E5f ) (--4 1) ) + E4 2) (--? ) )' ( 54 ) 
i=i j=i 

The resulting model is illustrated for the non polar case {Sj: — 0, and A"' s — Al' s V i) in Fig. O the parameters of 
the model are specified in the figure caption. An isolated atomic-like charge density is shown in the inset. The sizes 
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of the layers of the two materials that constitute the multilayer can be tuned by changing the number of building 
blocks in the left, N%, or in the right, N 2 . 

1.4 
1.2 
1.0 
0.8 
- 0.6 
0.4 
0.2 
0.0 
-0.2 

1 1 1 1 1 2 2 2 2 2 

FIG. 5: Illustration of the toy model used to simulate microscopic charge-densities in a multilayered material. A non-polar 
case is represented (6^ — and A"' s — A^' s V i). The remaining parameters have been chosen as A"' 1 = A e<1 — 1, a ny i — 0.5, 
cr e ,i = 2.0, A n ' 2 = A e ' 2 = 2, cr„,2 = 0.7, a e , 2 = 4.5, 5 W = <5 (2) = for all the "atoms", a w = 6, and a (2) = 8. Within each 
period of the multilayer, there are 10 atomic layers of the left material and 10 of the right. Only a small portion, centered at 
the interface, of the microscopic charge density of the supercell is shown here. Inset: charge density of an individual "atom" 
on the left. Atomic units are used. 



p (z) is formed from p (z) by convoluting it with the filter function. Such convolutions are most conveniently formed 
by fast Fourier transforms, which require the discretization of space into a uniformly-spaced grid of points. 

VII. RESULTS: TOY MODEL, NON POLAR CASE. 

In Fig. EJa) we illustrate the specific toy model we shall analyze in detail for the unpolarized case. The parameters 
of the model are A 71 ' 1 = A 6 ' 1 = 1, a nA = 0.5, a eA = 2.0, A n < 2 = A e ' 2 = 2, a n . 2 = 0.7, cr e , 2 = 4.5, 5« = = for 
all the "atoms", = 6, and a^> — 8. Atomic units are used throughout the paper. First, we consider a reference 
density p (z) defined as in Eq. ([6]), locating the interface at Zi„ t half way between the rightmost atomic layer of 
material 1 and the leftmost atomic layer of material 2. A~p u (z), defined as in Eq. ((7]), is plotted in Fig. [6]Jb). Next, 
we construct ~p u (z) from ~p u (z) directly using a nanosmoothing function /(z — z ) defined as in Eq. (|15p . in which 

11 = aS 1 ' = 6, and l 2 = = 8. The smoothed charge density ~p u (z) is shown in Fig. [6fc). As we have already 
discussed in Sec. IV CI the regions R s within which ~p u vanishes lie within the regions R s within which A~p u vanishes, 
because smoothing ~p u within R s brings into ~p u values of ~p u for z outside R s . Due to the symmetry of the microscopic 
charge density at the center of each layer, the interfacial charge density between material 1 and material 2 is the 
mirror image of the interfacial charge density between material 2 and material 1 [see the center and the edges of the 
Fig. Mb) and Fig. M.c)}- 

A closer look at the interface region is shown in Fig. \7\ The discontinuity of A~p u at the interface plane is clearly 
observed. Also, we can see how A~p u displays large fluctuations at the atomic scale in the neighborhood of the 
interface due to the interface-induced relaxations of the atomic layers. The positions of the atoms in the layers close 
to the interface do not coincide with the positions of the atoms after the cleaving of the bulk to define Po ( see Sec. 
IIV[) . Therefore, in the computation of A~p u we are subtracting charge densities centered on different positions. The 
nanosmoothing procedure eliminates not only the contributions associated with ~p in the bulk regions R\ and i?2, but 
it filters out the oscillations due to the interface induced relaxations while producing a continuous charge distribution. 
Nanosmoothing ~p u is clearly superior to forming A~p u . 

In Fig. [51 we show ~p u (z) for three different smoothing functions, defined as in Eq. (HD with h = aW = 6, and 

1 2 = a (2) = 8 ; h = 2aW = 12, and l 2 = 2a^ = 16 ; and l x = 3aW = 18, and l 2 = 3a^ = 24, respectively. The 
widths of the regions where ~p u (z) vanishes (R x and R 2 ) depend on the width of the filter function. As a general rule, 
the more extended the filter function, the narrower the R regions. Although the profiles of the charge density differ 
significantly, the three charge distributions have the same net charge density Q u and dipole moment density p u (see 
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FIG. 6: (a) Microscopic charge density ~p u (z) within a unit supercell for an unpolarized system. The widths of the left (1) 
and right (2) materials have both been set to 10 atomic layers, with interatomic separations of a' 1 ' = 6 and = 8. The rest 
of the parameters of the toy model are as in the caption of Fig. [5] (b) Interface-induced charge density Ap u (2) defined by 
subtracting the charge density of a reference system, denned as in Eq. (|6]). The interface plane Zint is positioned at the middle 
of the separation between the last atomic plane on the left and the first atomic plane on the right, (c) Nanosmoothed charge 
density ~p u (z) of the microscopic charge density. The nanosmoothing function is defined as in Eq. (115[l with h = tr 1 ' = 6, and 
I2 — a' 2 ' = 8, single interplanar distances along z for each material. Atomic units are used. 



Table [TJ . This striking sensitivity of the nanosmoothed charge density to the smoothing function impedes detailed 
physical interpretation of its features. 

Even though we are dealing with an unpolarized interface made of the juxtaposition of neutral and non-polar atoms, 
the interface charge density Q u , Eq. ()4a|) . vanishes if and only if the integration limits z\ and 22 are taken midway 
between atoms inside each material layer, meeting the requirements for symmetry described in Sec. IV El (see Table 
UJ. AQ U , Eq. (|8a[) . on the other hand, does not vanish since, with our criterion for locating the interface plane Zi n t, 
we do not conform to the requirement that Zi nt — z\ contains an integer number of unit cells of the left material, and 
Z2 — Zi n t contains an integer number of unit cells of the right material. Consequently, the integral of the reference 
charge density Po between z\ and zi carries a net charge. On the other hand, as it should, Q u vanishes and is 
independent of the integration limits provided that Z\ lies in R 1 and z% lies in R 2 ■ 

We have also calculated the corresponding interface dipole densities p u , Eq. (|4b|) ; Ap u , Eq. (|8b[) ; and p u , Eq. (|36b[) . 
with the results displayed in Table [U By definition, p u and Ap u are independent of the nanosmoothing procedure. 
As long as l\ and I2 equal an integer number of lattice constants of the bulk unit cell of the material along z, p u 
is insensitive to the shape and range of the smoothing function. Note that, as expected after the discussion in Sec. 
IV El the nanosmoothed dipole density p u differs from both Ap u , computed from Ap u (z), and p u computed from 
~p u (z). The latter two also differs because ~p (z) defined in Eq. (|6|) generates an extra dipole-moment density when 

integrated within our integration limits. Thus the electrostatic potential shifts AV U and AV U differ correspondingly. 

The shift AV U is the physically meaningful one because what enters in the band offsets are the local averages of the 
electrostatic potentials in each material, which are independent of the location of z\ and 22 within R Y and R 2 entering 

in Eq. (|27jl defining AV U - The analogous relation AV U , 

AV U = V u (z 2 ) - V u (zi) , (55) 

shows that AV is the difference of the two potentials at specific points within regions in which V (z) varies rapidly. 
In Fig. [5] we show V (z) and V (z) and indicate the positions z\ and Z2- It is clear that AV is about 20% larger 
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FIG. 7: Blow up of the interface-induced charge densities obtained from the subtraction of the reference bulk density and from 
nanosmoothing. Their difference is clearly shown. 
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FIG. 8: Nanosmoothed charge density p u (z) obtained by nanosmoothing the planar average charge density p u (z) shown in 
Fig. EJa) with three different filter functions. They are defined following Eq. (|15[1 with l\ — a ( 1} = 6, and h = a (2) = 8 (solid 
line); h = 2a (1) = 12, and h = 2a (2) = 16 (dashed line); and h = 3a (1) = 18, and h = 3a (2) = 24 (dotted line), respectively. 



than AV, explaining the relation between p u and p u in Table U the difference arising from the way V is sampled by 
nanosmoothing and by the selection of Z\ and z 2 ■ 

In Table U we also show two cases where the nanosmoothing does affect the value p u . First, the value of p u 
departs from the correct value 2.157 when the range of the smoothing function, L in Eqs. (113a|) and (Il3bl) . is of the 
same magnitude as the width of the layer of one of the materials, as for l\ = = 30 and 1% = 5a' 2 ' = 40, so 

L = h + 12 = 70, slightly larger than the width of material 1, made of 10 layers with an interlayer distance = 6. 
Under these circumstances, it is not possible to define regions R x and R 2 where ~p u (z) vanishes, impeding a proper 
location of the integration limits z\ and zi. Second, the same occurs when the range of every filter function entering 
Eq. IT5|) does not equal an integer number of lattice constants of the bulk unit cell along z, as for l\ — = 3, and 

(2) 

a case, the charge density after nanosmoothing, p u (z), still shows large and rapid oscillations. 
In Table [TT] we report the amount of charge transferred from one material to the other, the screening length and the 
interface position computed with the method summarized in Sec. IVDI for the nanosmoothed charge densities plotted 
in Fig. [5] The negative cumulative charges are displayed in Fig. [TUJ As expected the value of q is sensitive to the 
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TABLE I: Interfacial charges and dipole moments of the microscopic interface-like charge density shown in Fig. |6ja). The 
interface charge Q u and the dipole-moment p u densities are denned respectively in Eqs. (|4a|l and (|4b|l and the integration 
limits 21 and 22 taken midway between the atoms at the center of each material layer. When a reference charge density of an 
ideal interface is subtracted, the corresponding charge AQ U and dipole moment Ap u axe defined as in Eqs. (|8a[l and (I8b[l with 
the same integration limits as before. Since the nanosmoothed charge density does not play any role in the definition of Q u , 

Pu, AQ U , and Ap u , these magnitudes are insensitive to nanosmoothing. Q u and p u are defined in Eqs. (|36a[) and (|36b[l . l\ and 
h are the lengths of the square- wave functions entering in the definition of the smoothing function, Eq. (|15[1 . 



h 




h 




Qu AQu 


Qu 


Pu 


Ap u 


Pu 


= 


6 


a™ = 


8 


-0.113 





1.749 


-4.821 


2.157 


2a« = 


12 


2a^ = 


16 











2.157 


3a« = 


18 


3a' 2 ) = 


24 











2.157 


4a« = 


24 


4a< 2 > = 


32 











2.157 


5a« = 


32 


5a< 2 > = 


40 











2.121 


2 


3 


J 2 '' 

2 


4 











1.871 



30 



20 - 



- 10 - 



a 



- 



-10 



-20 



1 1 1 1 1 

— Vu (z) 


1 1 1 1 1 


1 






1 


1 






1 


"-- Vu{z) 
























'V 


V 


V 


\l 


f v 


V 




\- 


\ 


I 

AVu I 


















\ 
l\ 


i 
i 

j 










AVu 






\\ J\ J\ l\ l\ 






1 






1 


1 






i 



11111 



1111122 2 2 2i2 2 2 22 



Zl Z 2 



FIG. 9: Planar averaged electrostatic potential computed by solving the one dimensional Poisson equation, Eq. ((2}, with the 
charge density p„ (z) shown in Fig. [6ja), V u {z), and by solving Eq. (|19p with the charge density p u (z) shown in Fig. EJc), 
V u (z). The integration limits z\ in R x and 22 in R 2 are indicated by dotted lines. The electrostatic potential shifts AV U 

[Eq. (|55[) ] and AV U [Eq. (|27|l ] defined as the difference of the corresponding potentials at the points z\ and 22 are also shown. 
Atomic units are used. 



smoothing function, decreasing with the increasing of its range. Accordingly, the screening length is also sensitive 
to the range of the smoothing function L, taking a value that is roughly the greater of l\ and l-i- Nevertheless, the 
position of the interface seems to be insensitive to the filtering function. 

The results displayed in Fig. [8] for ~p(z) demonstrate that the larger the width of the smoothing function, the more 
complex the spatial dependence of ~p{z) and the less it resembles a simple interface charge density. The introduction 
of multiple extrema is artificial and caused by excessive transfer of charge back across the interface by smoothing. 
This back transfer of charge is manifested clearly in the values of q in Table [TTJ Accordingly the smallest allowable 
widths l s — a Sj should be used for nanosmoothing. The resulting values of q and A are the best estimates which can 
be extracted by nanosmoothing. 

Now, as a final test of the insensitivity of the dipole-moment density and potential shift to the shape of the smoothing 
function, we shall simplify the toy model to one with a common atomic separation in both materials. This opens 
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FIG. 10: Negative cumulative charges defined in Eq. (|30a|l for the nanosmoothed charge densities shown in Fig. [8] We can 
estimate the magnitude of the charge transferred per unit area of the interface q from the supremum of this cumulative charge 
and the interface position from the position where this supremum appears. Numerical values are given in Table ITT1 



TABLE II: Magnitude of the transferred charge density, screening length and position of the interface computed from the 
nanosmoothed charge densities plotted in Fig. [8] h and I2 are the lengths of the square- wave functions entering in the definition 
of the smoothing function, Eq. (|15[l . Since the interface is unpolarized, the charge at each side of the interface must be equal 
in magnitude and oposite in sign. 



h 


I2 Zi n t 


q 


A 


6 


8 58.083 


0.253 


8.517 


12 


16 57.709 


0.143 
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three simple options for the filter function. We can set / equal to a single lu when the filter function adopts a square 
shape [see Fig. E£a)], we can set / equal to the convolution of two w's with the same I, a triangle function, 



/(*-*) 



(56) 



or we can set / equal to the convolution of two oj's with different I, a trapezoidal function, as in Eq. (|17p and in Fig. 
EJb). Comparing the results of the smoothing with the square wave, the triangle, and the trapezoidal function for 
various values of I allows us to establish which physical properties are insensitive both to the range and to the shape 
of the smoothing function. 

Again, as soon as l\ and/or I2 equal an integer number of lattice constants of the bulk unit cell of the relevant 
material along 2, p u is insensitive to the shape (triangular, square or trapezoidal) and range of the smoothing function. 



VIII. POLARIZED SYSTEMS. 



Up to now we have dealt with materials unpolarized except at the interfaces. In such systems it was possible to 
identify bulk-like regions in the middle of each layer where the microscopic charge density is unaffected by the presence 
of interfaces and returns to the bulk microscopic charge density, which is centrosymmctric within a unit cell. We now 
turn our attention to polarized interfaces, where at least one of the materials has a non-vanishing polarization P. 

Basic electrostatic arguments that can be found in any textboo k 21 ' 32 show that a non uniform polarization in a 
dielectric generates a volume charge density, the polarization charge p po i (r), whose value at any point of space is 
given by 
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TABLE III: Interfacial charge and dipole moment densities for a microscopic interface-like charge density obtained with the 
following parameters: A"' 1 = A"' 1 = 1, cr n>1 = 0.5, cr M = 2.0, A"' 2 = A e ' 2 = 2, cr„, 2 = 0.7, cr e , 2 = 4.5, 5 W = S (2) = for all 
the "atoms", a (1) = 10, and a (2> = 10. The meaning of the symbols and integration limits is as in Table U 
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p P oi(f) = -V-P(f). (57) 

Even in the case of uniformly polarized material, where P is constant and therefore its divergence vanishes inside the 
dielectric, the discontinuity of the polarization at the surface or interface gives rise to a net surface or interface charge 
density a po i given by the familiar form 



Vpoi = P • n, (58) 

where h is a unit vector normal to the surface or interface pointing outwards. 

In the case of an interface, the other material at the interface responds to this perturbation in order to minimize the 
electrostatic energy cost associated with the build up of the polarization charge V • P at the interfaces. If the second 
material is a dielectric, a uniform polarization is induced within it as wel h 33 i 34 If the second material is a metal, 
a screening charge is induced that spreads over a finite distance (the screening length) in the electrode, producing 
additional dipole layers at each polar dielectric/metal interface. If the screening of the polarization charge is not 
perfect, a residual depolarizing field appears inside the dielectric 

We now face the same question as for the unpolarized case: how to extract from the results of the first-principles 
calculations values of the physical quantities of interest (polarization and screening charge densities, screening lengths, 
depolarizing fields, etc). 

In parallel to what was done for the unpolarized case, we can define the net charge density Q p and dipole moment 
density p p associated with the microscopic charge density of the polar interface ~p p as 



Q p = dz p p (z) , (59a) 

Jz ± 

i>Z2 

p p — dz zp p (z) . (59b) 

J Z\ 

The computation of these two quantities present the same difficulties as for the unpolarized case, summarized in Sec. 

inn 

The approach of determining a reference interface charge density 7j (z) as in Eq. ([6]) and defining an interface- 
induced deformation of the charge density A~p p (z) by subtracting ~p (z) from the microscopic charge density, 



Ap p (z) = p p (z) - p (z) , (60) 

discussed in Sec. II VI for the unpolarized case, still allows for the identification of two regions R\ and i?2 where A~p p (z) 
vanishes [see Fig. ITTT b)]. Therefore we can define both the interface charge and dipolar densities invariant with 
respect the position of the integration limits z\ and Z2, as long as z\ lies in R\ and z-i lies in i?2- 

However, on top of all the drawbacks to that approach presented in Sec. IIVI there is an additional pitfall now: 
an extra dependence on the position of the integration limits Z\ and z 2 appears in the computation of the interface 
dipole moment from the definition of the reference interface position Zi nt . Indeed, for at least one of the materials the 
dipole moment density does not vanish within the bulk unit cell, 
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dz zp ( Q s > 0) . (61) 

Therefore, the dipole moment density Ap p associated with A~p p (z) does not equal p p defined in Eq. (|59b() . not even 
in the case where Zi nt — Z\ contains an integer number M\ of unit cells of the left material and z 2 — z int contains an 
integer number M 2 of unit cells of the right material, 

Ap p = dz zAp p (z) = dz z [p p (z) - p (z)] 

J Z\ J Z:\ 

= I dz zp p (z) — I dz z~pQ (z) 

= p p -M 1 p£ ) -M 2 pf\ (62) 

Nevertheless, the nanosmoothing procedure developed in Sec. IV Al remains a useful tool, since the methodology for 
constructing the nanosmoothed charge density, Eq. (fT2|) . is independent on whether the microscopic charge density 
is polarized or not. Both the nanosmoothed Poisson equation [Eq. (Tl9|) ] and the analysis of the insensitivity of the 
dipole moment density to the smoothing function (Sec. IV CP still hold under the same conditions for the smoothing 
function as for the unpolarized case, because the fundamental starting point for the derivation, the microscopic 
Poisson equation and the existence of regions with a negligible nanosmoothed charge R l and R 2 , are insensitive to 
the polarization state of the microscopic charge density. 




IX. RESULTS: TOYMODEL, POLARIZED. 

In Fig. [TTT a) we illustrate the first specific toy model we shall analyze in detail for the polarized case. The 
parameters of the model are the same as in the caption of Fig. O with the exception of the parameter 8[ , that 

now takes a non-zero value 5^' — 1.0 where i runs over all the "atoms" of material 1. This toy model simulates a 
multilayer in which a net polarization has been induced in material 1 by displacing the "electronic clouds" rigidly 1.0 
length units towards the right. The neutral but polarized atom of this toy model can be thought of as representing 
the overall neutral atomic planes of a ferroelectric within which a dipole-moment density is generated by buckling 
or puckering. In principle, the charge density of material 2 would be modified as a response to the presence of a 
polarization in material 1. This polarization- induced response, which should be computed self-consistently, is not 
considered in the present simple toy model. The displacement of the negative charges translates into the asymmetry 
of the charge density inside material 1 represented in Fig. Illf a). The interface-induced deformation of the charge 
density A~p p (z), defined as in Eq. (|6"01) , is shown in Fig. [TUb) where, as for the unpolarized case, the position of the 
interface Zi nt has been located half way between the rightmost left atomic layer and the leftmost right atomic layer. 
The nanosmoothed charge density ~p p (z) is displayed in Fig. Illf c). obtained with the same smoothing function as the 
the one used in the construction of ~p u (z) in Fig. E^c). In the last two panels, the two regions inside materials 1 and 
2 where A~p p (z) and ~p p (z) vanish are clearly shown. 

Some of the features already discussed in Sec. I VIII for the unpolarized interface still remain valid for the polarized 
case. In particular: (i) the regions R a within which ~p p vanishes lie within R s within which Ap p vanishes, because 
smoothing ~p p within R s brings into ~p p values of ~p p for z outside R s ; (ii) the discontinuity of A~p p at the interface 
and the large fluctuations at the atomic scale in the neighborhood of the interface; (Hi) the smooth and continuous 
profile [note the change of scale in the charge density in Fig. [TlTb') and Fig. HHc)] of ~p p , highlighting its advantages 
with respect Ap p . However, some other issues are new. Now, as a consequence of the asymmetry of ~p p (z), the charge 
density at the interface regions, defined as the ranges where Ap p (z) and ~p p (z) differ significantly from zero, are 
different in the adjacent interfaces contained in our simulated supercell [see the center and the edge of Fig. HHh) and 
ITTT c)]. They are not simply reflections of one another as in Fig. [6jb) and[6]Jc). 

In Table HVl we report the values of the interface charges computed from p p (z) (Q p ), Ap p (z) (AQ P ), and p p (z) 

(Q p )- The integration limits z\ and z 2 are taken midway between the atoms at the center of each material layer. In 
contrast to the unpolarized interface, Q p does not vanish. The surplus of negative charge is due to the charge density 
entering region z 2 — Zi nt from the left from the region between Zi n t — z\ which is not compensated by the departure 
of any charge. Since ~p p (z) is not even about z%, Q p ^ Q p . As happened in the unpolarized case, as long as l\ and l 2 
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FIG. 11: (a) Microscopic charge density p p (z) within a unit supercell for a polarized system, with the polarization induced by 
a displacement of the electronic charge in one of the materials. The widths of the layers of both material and the rest of the 
parameters of the toy model are as in the caption of Fig. \5\ with the exception of 5^ = 1.0, where i runs over all the "atoms" 
of material 1. (b) Interface-induced deformation of the charge density Ap p (2) defined as in Eq. (|60p . As before, the interface 
plane Zmt is positioned at the middle of the separation between the last atomic plane on the left and the first atomic plane on 
the right, (c) Nanosmoothed charge density p p (z) of the microscopic charge density. The nanosmoothing function is defined as 

in Eq. (|15[) with h = a' 1 ^ = 6, and h = a' 2 ' = 8, single interplanar distances along z for each material. Atomic units are used. 



equal an integer number of lattice constants of the bulk unit cell of the material along z, Q p and p p are insensitive to 
the shape and range of the smoothing function. Only when the range of every filter function entering in Eq. (|15j) does 
not equal an integer number of lattice constants of the bulk unit cell, or when the range of the smoothing function L 

is of the same order as the width of the layer of one of the materials, do Q p and p p differ from the correct value. 

Subtracting from ~p p (z) the related p u (z), a microscopic charge density constructed by nanosmoothing with filter 
functions defined with the same parameters with the exception of the displacement of the electronic charge, we obtain 
the profile of the nanosmoothed polarization-induced charge density, Fig. 1121 whose integral between Z\ and Z2 gives 

<jpoi as the accumulation of charge at the interface. Note that since the resulting Q u vanishes, Q p is identical to a po i. 
Although, as listed in Table IIV1 the amount of charge accumulated at the interface is independent of the shape and 
range of the filter function used (provided that l\ and I2 are integer numbers of lattice constants of the bulk unit cell 
of the material along z), the profile of the polarization-induced charge density is not. The wider the nanosmoothing 
function, the larger the range of the polarization charge density. Again, this sensitivity of the shape of the smoothed 
charge density to the nanosmoothing function spoils a direct physical interpretation of its features. The smallest 
values of l\ and 1%, a single lattice constant of each material, yields, as for the unpolarized case, the simplest and 
physically most relevant profile for ~p p . 

Solving the nanosmoothed Poisson equation, Eq. (flT?)) . for the nanosmoothed charge density of Fig. [TTT al we obtain 

the nanosmoothed potential V p (z) displayed in Fig. [13] This electrostatic potential has two distinct features: (i) a 
jump at each interface due to the interface dipole moment already present in the unpolarized case (see Fig. [9] and 
TableJT]); (ii) a field generated by the polarization charge a po i at the interface. We can isolate this field by subtracting 
the nanosmoothed potential for the polarized (Fig. [T3"]) and unpolarized (Fig. \§§ systems. The result is shown in Fig. 
IT4l The electric field can be computed from the slope of the nanosmoothed potential. 

The magnitude of the electric field is related to the magnitude of the polarization charge u po \. The periodic 

boundary conditions enforces that £^Nia^, k + ^ 2 'JV2a^L = 0, where N% and N2 are the number of unit cells of 
materials 1 and 2 stacked to build the superlattice (see Sec. IVII) . Since each interface carries a charge of magnitude 
<7 po i, the electric field changes its value at the interface as E ^ — S^- 1 ' = 47r<7 po ;, yielding to values for the fields 
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FIG. 12: Difference between the nanosmoothed charge density for the polarized interface and an unpolarized interface. The 
parameters of the microscopic charge density as are in Figs. [TT]and[5l respectively. Three different filter functions are used for 
the nanosmoothing, with the same parameters as in Fig. [8] 
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FIG. 13: Planar averaged electrostatic potential computed by solving the one dimensional Poisson equation, Eq. P). with the 
charge density p p (z) shown in Fig. Illf a'l. V p (z), and by solving Eq. (|19|l with the charge density ~p p (z) shown in Fig. Illf c). 

V p (z). The nanosmoothing function is denned as in Eq. CEI with h = a (1) = 6, and h = a {2) = 8, single interplanar distances 
along z for each material. Atomic units are used. 



£(i) = jff PoiN^a and g (2) = ^ P oiN ia _ Taking a vol from Table|lVl we infer a value of 8& = 1.196, and 

£^ — —0.897, in excellent agreement with the slopes computed in Fig. [14] 

From the knowledge of the change of the macroscopic electric field across the interface we can determine the 
difference in the zero-field polarization of the two materials. Assuming that in our toy model the dielectric constant 
of the two materials are the same and equal to 1, then 



Since £^ - = Aira poll 



-47T (v^-V^) =e(eW-£W) 



(63) 
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FIG. 14: Difference of the nanosmoothed potential for the polarized (dashed line in Fig. I13[) and unpolarized (dashed line in 
Fig. [9j) cases. The magnitude of the electric field £ generated by a periodic array of polarization charge distributions tj po i can 
be computed from the slope of the profile of the potential within material 1 and material 2. Atomic units are used. 



TABLE IV: Interfacial charges and dipole moment densities of the microscopic interface-like charge density for the polar 
interface shown in Fig. Illf a) . The meaning of the different symbols is as in the caption of Table [I] with the subscript u 
replaced by p for the polar case. 
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-(pM-V™) =a pol . (64) 
But our material 2 is unpolarized, so = 0, and we arrive to the conclusion that 

V {1) =a po i. (65) 

This conclusion can be checked analytically in our particular toy model, where the charge density is the juxtaposition 
of "atomic-like " charge densities. Integrating the first-moment of the microscopic charge density for the slab of 
material 1 used to build the superlattice, defined in Eq. I53|) . we arrive to the conclusion that = — ttj~, whose 

a bulk 

numerical value in our numerical example equals the polarization charge shown in Table lIVl A generalization for the 
dielectrically mismatched interface can be found in Sec. III. F of Ref. [35|. 

A second toy model we shall analyze for the polarized case is shown in Fig. [TH It is made of a one dimensional 
chain in which the unit cell of material 1 has two atoms per unit cell. Inside the unit cell, one of the atoms is charged 
negatively and the other positively so that the overall charge in the unit cell is neutral. Thus, besides the electronic 
polarization discussed in the previous case, polarization can be induced by ionic displacements. The parameters of the 
model are as in Fig. [Til with the exception of A e ^ = 1.5 for the odd "atoms" (net charge, -0.5 per site), A e ' 1 = 0.5 

for the even "atoms" (net charge, +0.5 per site), and 8^ = 0.0, where i runs over all the "atoms" of material 1 (no 
displacement of the electronic cloud is considered). 
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FIG. 15: (a) Microscopic charge density p p (z) within a unit supercell for a polarized system in which the unit cell of material 
1 has two atoms of opposite charge in the unit cell. The widths of the layers of both material and the rest of the parameters 
of the toy model are as in the caption of Fig. 1111 with the exception of A"' 1 — 1.5 (net charge, -0.5 per site) for the odd atoms 
and A e ' x = 0.5 (net charge, +0.5 per site) for the even atoms of material 1, and = 0.0, where i runs over all the "atoms" 
of material 1. (b) Interface-induced deformation of the charge density Ap p (2) defined as in Eq. (|60[) . As before, the interface 
plane Zi nt is positioned at the middle of the separation between the last atomic plane on the left and the first atomic plane on 
the right, (c) Nanosmoothed charge density ~p v {z) of the microscopic charge density. The nanosmoothing function is defined as 
in Eq. (|15fl with l± = = 6, and I2 = e^ 2 ' = 8, single interplanar distances along z for each material. Atomic units are used. 



The main conclusions that can be drawn from Fig. [15] are the same as in Fig. 1111 As we made for the electronic 
polarized case, we can infer the value of the bulk zero field polarization from the polarization charge at the interface, 
that in this case amounts to +0.25. 



X. SUMMARY AND DISCUSSION. 



The problem of how to extract from the immense detail provided by the first-principles calculations, with resolution 
at the atomic scale, reliable values of physical quantities of interest which enter into nanoscale electrostatic analysis 
has been reviewed. 

This problem is particularly challenging in the case of interfaces between quite different materials, since all the 
relevant magnitudes (interface charge and dipole densities, screening lengths, etc.) are overwhelmed by the large and 
rapid oscillations of the microscopic charge density. 

The different procedures to filter out the periodic oscillations of the microscopic quantities, which typically follow 
the underlying atomic structure, preserving only those features that change in the vicinity of a surface or interface 
are critically analyzed, and the criteria under which they accurately extract the quantities of interest are discussed. 

The approach of defining a reference interface charge density from the bulk unperturbed charge densities of each 
material is spoilt by the fact that the interface position is undetermined. The profile of the charge density at the 
interface is discontinuous and displays large fluctuations in the neighborhood of the interface due to the interface- 
induced atomic relaxations. Moreover, both the interface charge and dipole densities (and so the potential drop at 
the interface) are not unique in this approach, since they depend critically on the position of the interface. 

A clearly superior method consists in nanosmoothing the microscopic charge density by taking its convolution with 
a filter function. This procedure eliminates the contributions coming from the bulk together with the oscillations due 
to the interface-induced relaxations while producing a continuous charge distribution. 

In this work, we prove rigorously that the interface charge and dipolar densities are independent of the nanosmooth- 
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ing function used provided the following conditions are met: (i) it is positive definite, even and normalized within the 
support region of space within which it is defined [Eqs. (|13a,[l - (|13d[) ] ; (ii) it must be sufficiently smooth itself for the 
Poisson equation to remain invariant after nanosmoothing. By "sufficiently smooth" we understand that the second 
derivative must exist, and both the function and the first derivative must vanish at the end points. [Eqs. (|18a[) - (|18c[) ]; 
(Hi) in practice the convolution is done using fast Fourier transforms. Therefore, although the nanosmothing function 
might violate some of the previous conditions, the family of the finite Fourier series involved in the convolutions is 
a distribution which converges to the smoothing function in the limit, meeting all the requirements along the way. 
(iv) if square-wave filter functions (or convolutions of them) are chosen as filtering functions, the width of each filter 
function l\ and I2 must be integral numbers of the lattice constant of each material in order to have charge and dipole 
moment densities (and therefore potential shifts across the interface) insensitive to the smoothing function; (v) the 
smallest acceptable value for li and l 2 should be used in order to avoid the presence of multiple maxima and minima in 
the smoothed charge density; (vi) the total width of the filter function L should be chosen on the scale of the unit cell 
length or larger, but significantly smaller than the regions R s of each material where the microscopic charge density 
is unaffected by the presence of the interface; (vii) after nanosmoothing, the charge density displays two regions R s , 
within the R s , where the smoothed charge density ~p vanishes. The integration limits z\ and z 2 used to compute the 
interface charge and dipole densities must lie in these regions; (viii) once the integration limits are chosen, the charge 
density at the interface computed from the microscopic and the nanosmoothed charge density are equal if and only 
if the microscopic charge density is symmetric for \z — z\^\ < L. The interface dipole density computed from the 
microscopic and the nanosmoothed charge density are never equal; (ix) nanosmoothing is only valid for computing 
charge and dipole moment densities. Nothing can be said about the shape of the nanosmoothed charge density at 
the interface since it depends critically on the filter function; (x) by using the smallest acceptable values of l\ and l 2 , 
generally one unit cell length, reasonable estimates of the density of the charge transferred across the interface and 
of dipole-layer width can be made. 

The nanosmoothing procedure is a powerful technique that allows us to extract the information relevant for comput- 
ing the change in the average potentials and charge densities from one side of the interface to the other, opening the 
door to the calculations of band offsets^^ ' 17 i 18 polarization 29 and dielectric permittivity^ profiles, effective charges^ 
and force constants^ in semiconductor-semiconductor interfaces, and depolarizing electric fields and screening lengths 
in real ferroelectric capacitorsi^ The present work sets the basis for a rigorous selection of the minimum size of the 
supercell required to obtain accurate and smoothing-independent results. 
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APPENDIX A: TRANSVERSE AVERAGING OF THE POISSON EQUATION. 

The microscopic Poisson equation is written as 

V 2 V(r) = -4np(r). (Al) 
Integrating both sides of Eq. (IA1|) in the (x, y) plane and dividing by the surface S of the interface unit cell, we get 

~S S S V2y ^ dxdy = -4w^J J p(r) dx dy 

= -Aifp{z). (A2) 
We now write out the integral of the left-hand-side of Eq. (| A2[) explicitly, 



1 f f d 2 V(r) , , 1 f fd 2 V(r) , , 1 f f d 2 V (r) , , 

' ' ' dxdy+- ; ' dxdy+- / ; ' dx dy. (A3) 



SjJsdx 2 " S J Js dy 2 a S J J s dz 2 
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The first two integrals on the right-hand-side can be performed trivially, 



d 2 V (f) , , If dV (r) 

^^ dxdy = s ~^x~ 



dy = 0, 



(A4) 



boundary x 



where boundary x refers to the two intersections of the boundary of the unit cell with the x axis. Since the potential 
and all its derivatives are periodic in the plane, the previous integral vanishes. The same holds for the second integral 
in the right-hand-side of Eq. (|A3[) . 

Regarding the third integral in Eq. (|A3[) . we can take the second derivative with respect z out of the integral, 



d 2 V(r) , , d 2 (1 
2~ dx d V = 7T2 



dz 2 \S j Js 
d 2 V(z) _ d 2 V(z) 



V (r) dx dy 



dz 2 dz 2 

Gathering together the results of Eqs. (|A2j) through (|A5p we can conclude that 



(A5) 



V 2 V (r) = V 2 V (z) 



d 2 V(z) 
dz 2 



-4np(z) ; 



(A6) 



that is, the transverse average of the Poisson equation yields a one-dimensional Poisson equation for the transverse 
average of the potential. 

APPENDIX B: NANOSMOOTHING OF THE POISSON EQUATION. 

Applying the nanosmoothing procedure to the electrostatic potential V (r) results in 



V(z)= / dz f(z-z')V(z), 
where V (z) is the planar average of the of V (r) in planes parallel to the interface, 



(Bl) 



v{z) = s J V ^ dxdy ' 



(B2) 



and f(z — z ) vanishes for z > z + L and z < z — L. Taking the second derivative of the nanosmoothed potential 
and requiring that d ^a" z - exists yields 



V'V(z) = 



dz 2 

Z+L , , d 2 
L 

z+L j2 

dz ( 
z-L dz- 

L , , d 2 f{z-z ) —. ,. 
dz — V(z ) 



f(z-z)V(z) 
d2 ^- Z '±V(z') 



-L 



df(z - z ) 
dz' 



dz 



no 



,2 



z+L 



z+L 

dz 

-L dz 



df(z-z') dV(z') 
dz' 



(B3) 



Requiring that 
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df(z) 
dz 



df(z) 
dz 



= 0, 



(B4) 



makes the first term in the right hand side of Eq. (|B3j) vanish. Integrating the second term in the right hand side of 
that equation yields 



2+1 ' dz' ~ dV ( Z '^ 

-L 



dz' 



dz' 



dV{z 



dz 



z+L 



+ I dz f(z- z ) — r7T - 



— L 



dz' 



(B5) 



so that we arrive at the conclusion that 



V 2 V(z) 



dz 2 



dz f(z — z ) 

dz' 



d 2 V(z ) d 2 V(z) 



dz 2 



(B6) 



APPENDIX C: ELECTROSTATIC POTENTIAL SHIFT AND INTERFACE DIPOLE DENSITY. 



The formal solution of the Poisson equation, Eq. (fl~9|) . is 



V(z) = -2ir dz \z - z | p(z ) 



-2tt 



dz (z — z) p(z ) + / dz (z — z) p(z ) 



The electrostatic potential shift across the interface is, therefore, 



(CI) 



AV =V(z 2 )-V( Zl ) 
—An I dz z p(z ) 



' Zl 

2tt{ z 2 



dz p(z ) 



dz p(z 



dz p(z ) — dz p(z ) 



(C2) 



The location of the points at ±oo are thus far unspecified. They can both be chosen at an image of z\ or of z 2 
within the unit cell without loss of generality since all the relevant computational procedures require periodicity. In 
that case, all four integrals within the square brackets of Eq. (|C2[) vanish because of overall electrical neutrality and 
of the choice of Z\ and z 2 at the microscopic unit cell boundaries, with Eq. (|2H)) following. 
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